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Abstract 

We have developed a transferable nonorthogonal tight-binding total en- 
ergy model for germanium and use it to study small clusters. The cohesive 
energy, bulk modulus, elastic constants of bulk germanium can be described 
by this model to considerably good extent. The calculated bulk phase di- 
agram for germanium agrees well with LDA results. The geometries and 
binding energies found for small Ge^ clusters with n = 3 — 10 are very close 
to previous ah initio calculations and experiments. All these results suggest 
that this model can be further applied to the simulation of germanium cluster 
of larger size or with longer time scale, for which ah initio methods is much 

more computational expensive. 
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In the past decade, tight-binding molecular dynamics (TBMD) has evolved into a power- 
ful approach in the simulation of semiconductor materialsiHl. In the tight-binding scheme, 
although the system is still described in a quantum-mechanical manner, the computational 
cost is significantly reduced due to the parameterization of Hamiltonian matrix elements. 
In many cases of material simulations, it might offer a satisfactory compromise between 
empirical^ and first principleiJ^ methods for modeling the interatomic interaction. As an 
alternative of the accurate but costly ab initio molecular dynamics, TBMD can handle more 
complicated systems with acceptable accuracyEHi. 

For carbon and silicon, there are several well established orthogona and 
nonorthogona tight-binding models. Although the orthogonal models works well for 
various bulk systemsEHi, Menon found that the inclusion of the nonorthogonality of tight- 
binding basis is essential for describing the geometries and binding energies of small silicon 
clusterffl Compared to carbon and silicon, there is much fewer tight-binding models de- 
veloped for germanium. Recently, M. Menon has extended the nonorthogonal tight-binding 
(NTB) scheme to germanium and calculated the structures and cohesive energies of small 
Ge„ clustersEl. Although the cluster geometries obtained in Ref. [13] are generally consistent 
with ab initio results, the binding energies are overestimated. In this work, we perform an 
independent fitting of NTB parameters for germanium, which describes binding energies of 
germanium clusters better than that in Ref. [13]. This model is employed to study some bulk 
properties and considerably good results are obtained. 

In Menon's NTB mo delM, the total binding energy E/, of a system with A^^^ atoms can 
be written as a sum 

Eb = Eel + Ej-ep + NaEo (l) 

Eel is the electronic band energy, defined as the sum of one-electron energies for the 
occupied states: Eei = YlT'^^k- In Eci-(1)5 a constant energy correction term NaEo and a 
repulsive interaction E^ep are also included. 

On nonorthogonal basis set, the eigenvalues of system are determined from the secular 



equation: 

det\Hij - eSijl = 0. (2) 

Here the overlap matrix elements Sij are constructed in the spirit of extended Hiickel 
theory0, 

and the nonorthogonal Hamiltonian matrix elements by 

H^, = V.A^ + ^-Si] (4) 

where 

b2 = [b) 

is the nonorthogonality between two sp^ bonding orbitals and if is a environment dependent 
empirical paramete 

The Hij and Sij depend on the interatomic distance through the universal parameters 
Vij, which are calculated within Slater-Koster's scheme^. The scaling of the Slater-Koster 
parameters V\y^ is taken to be exponential with the interatomic distance r 

Vxy,{r) = V,y,{d,)e-''^'-'°^ (6) 

where do = 2.45A is the bond length for germanium crystal in the diamond structured. 

The repulsive energy Erep in Eq.(l) is given by the summation of pairwise potential 
function 

Erep = E E Xin,) = E E Xoe-^'^(--'^») (7) 

i j>i i j>i 

where Vij is the separation between atom i and j. 

In practice, we adopt the Slater-Koster hopping integrals V\A>((io) fitted from the band 
structure of bulk germaniumS. The on-site orbital energies e^, are taken from atomic 



calculationsEi. The only four adjustable parameters a, K, xo, Eq are fitted to reproduce the 
fundamental properties of germanium bulk and dimer. The input properties include: the 
experimental values of bulk interatomic distance 2.45 A0, dissociation energy 2.65 eV0 and 
vibrational frequency (286±5 cm~^)0 of Ge2 dimer, as well as theoretical bond length of 
Ge2 (2.375 A) from accurate quantum chemistry calculation at G2(MP2) levels. The fitted 
parameters are given in Table I. 
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TABLES 



Table I. Parameters in the NTB mode developed for germanium in this work. See text for 
detailed descriptions. 





ep 


^SScr 


^spcr 




^ppn 


-14.38 eV 


-6.36 eV 


-1.86 eV 


1.90 eV 


2.79 eV 


-0.93 eV 




do 


K 


Xo 




a 


Eq 


2.45 A 


1.42 


0.025 eV 




1.748 


0.79 eV 



We can first check the validity of the present NTB scheme by studying the fundamen- 
tal properties of germanium solid in diamond phase. The obtained cohesive energy 3.58 
eV/atom is very close to experimental value 3.85 eV/atom@. Furthermore, we have calcu- 
lated the bulk modulus B and elastic constants Cu, C12, C44 of germanium and compared 
with experimental valued! in Table II. Most of the bulk elastic properties such as B, Cu, 
C12 are well reproduced except that the C44 is overestimated by 0.35 Mbar in our model. 



Table II. Bulk modulus and elastic constants (in units of Mbar) of bulk germanium in 
diamond structure. NTB are the theoretical results from present NTB model; Exper. denote 
the experimental values taken from Ref.[22]. 





Cii 


C12 


C44 


B 


NTB 


1.125 


0.545 


1.019 


0.738 


Exper. 


1.288 


0.483 


0.671 


0.751 
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FIG. 1. Cohesive energies vs. relative volume for bulk germanium in simple cubic (sc), diamond 
and clathrate I phase from NTB (solid line) and LDA (open circle) calculationsil. 

By using the NTB scheme, we have also calculated the equation of states of germanium 
in different phases. In Fig.l, we present the zero-temperature phase diagram of the fee, sc, 
diamond and type I clathrate obtained from NTB model, along with recent LDA plane- 
wave pseudopotential calculationsil. It is worthy to noted that our NTB model is able to 
described the energy and atomic volume of clathrate phase. The energy of clathrate I is 0.06 
eV/ atom higher than that of diamond phase and its relative volume is about 15% larger than 
diamond phase. These results are consistent with the 0.08 eV energy difference and 15% 
volume change from LDA calculationii. The success in clathrate phase is important since the 
clathrate is also four-coordinated structure^i. On the other hand, it is natural to find that 
the agreement between LDA and NTB scheme become worse in the high-coordinated phases 
like fee since the present model is fitted from the diamond phase and dimer. However, the 
relative poor description of high coordinate phase will not influence the study on germanium 
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clusters since such high coordination (~ 12) does not exist in the geometries of germanium 
clusters. Considering its extreme simplicity and small number of adjustable parameters, the 
current NTB scheme gives a sufficient satisfactory overall description of bulk germanium 
properties. Therefore, one can expect that the model to give a reasonable description on 
the germanium clusters. 

In this paper, we determine the lowest energy structures of the Ge„ clusters with n = 
3 — 10 by using TBMD full relaxation. The ground state structures of Ge„ (n = 5 — 10) 
are presented in Fig. 2 and the geometrical parameters of small Ge^ (n = 3 — 7) clusters are 
compared with previous ab initio calculationslUHl in Table 111. In general, both the lowest 
energy structures and their characteristic bond length agree well with ab initio results. A 
brief description is given in the following. 




FIG. 2. Lowest energy structures of Ge^ (n = 5 — 10) clusters. 



Table III. Lowest energy geometries (with characteristic bond length parameters) of small 
Ge„ clusters. Tight-binding calculation (NTB) are compared with previous ab initio results 



such as: MRSDCiS, B3LYpi, LDAllil. The label of atom and bond for Ge„ are taken 
from Ref.29. 



n 


Sym. 


Bond 




Bond length (A) 












MRSDClii 


B3LYpil LDAil 


ldaS 


NTB 



3 C2v 1-2 3.084 3.070 3.20 2.91 2.71 
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D2h 



Oh 



5h 



1-3 
1-2 
1-3 
1-2 
1-3 
3-4 

1- 2 

2- 3 
1-2 
1-3 

3- 4 



2.320 
2.477 
2.622 
3.277 
2.456 
2.456 



2.312 
2.475 
2.619 
3.135 
2.476 
3.320 
2.553 
2.941 



2.26 
2.40 
2.53 
3.19 
2.39 
3.19 
2.47 
2.85 
2.65 
2.57 
2.59 



2.21 
2.35 
2.44 
3.10 
2.34 
3.10 
2.40 
2.78 
2.56 
2.49 
2.51 



2.38 
2.44 
2.57 
2.87 
2.44 
3.40 
2.47 
2.70 
2.83 
2.57 
2.53 
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The minimum energy structure found for Ge^ is an isosceles triangle (C2v) with bond 
length 2.38 A and apex angle 6 = 69.5°, in agreement with ab initio calculations (see Table 
III for comparison). The linear chain has higher total energy of about 0.95 eV. 

The ground state structure of germanium tetramer is a planar rhombus (D2h) with side 
length 2.44 A and minor diagonal length 2.57 A. This structure has been predicted as 
ground state in all ab initio calculationsi^&ii and the tight-binding bond length are close 
to ab initio results. 

For the Ges, the lowest energy configuration is obtained as a strongly compressed trigonal 
bipyramid (Dsh). The energy of structure is lower than the perfect trigonal bipyramid by 
0.62 eV and the planar edge capped rhombus by 0.35 eV. The trigonal bipyramid structure 
has been considered in all of the previous ab initio studieSffii. In those LDA based 
simulation without symmetry constraintilSS, the trigonal bipyramid is found to undergo 
severe compression and relax to the structure in Fig. 2. 

A distorted octahedron (D4h) is obtained for Gcg as lowest energy structure. This 
structure is found to be energetically degenerated with a edge-capped trigonal bipyramid 
{AE = 0.018 eV). This result agree well with recent B3LYP and Car-Parrinello calculation 

of Gejiii. 

In the case of Gey, we find a compressed pentagonal bipyramid with Dsh symmetry 
as ground state and energetically lower than the face capped octahedron by 0.63 eV. The 
pentagonal bipyramid structure has also been obtained from LDA based simulationSSE^. 



Table IV. Binding energy per atom Eh/n (eV) of Ge„ clusters obtained within the present 
NTB model, compared to experimental valueSol, ab initio results based on G2(MP2) leveS 
or LDA plane-wave pseudopentiafl'S, as well as nonorthogonal tight-bindingjii calculations. 



n 


Exper.Eiii 


G2(MP2)il 


ldaS 


LDAi 


NTB0 


NTB (present) 


2 


1.32 


1.25 


1.89 




1.31 


1.32 


3 


2.24 


2.02 


2.78 


2.66 


2.11 


2.06 
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4 2.60 2.49 3.32 3.19 2.66 2.62 

5 2.79 2.68 3.58 3.45 2.85 2.73 

6 2.98 - 3.76 3.63 3.05 2.95 

7 3.03 - 3.90 3.77 3.19 3.09 

8 3.04 - 3.82 3.69 3.17 3.05 

9 3.04 - 3.93 3.79 3.25 3.12 

10 3.13 - 4.04 3.91 3.32 3.17 
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An additional atom capped to pentagonal bipyramid of Gej yields the lowest energy 
structure for Geg. This structure is more stable over the bicapped octahedron by 0.08 eV. 
Both of these two structures are found for Ges in Car-Parrinello simulation, while bicapped 
octahedron is lower in energy by 0.03 eVii. 

A bicapped pentagonal bipyramid is found for Geg. It is more stable than a capped 
distorted square antiprism by 0.06 eV. The current ground state structure has been found in 
Car-Parrinello simulation for Geg^ but it is 0.08 eV higher than the capped square antiprism 
structure. 

For Geio, the tetracapped trigonal prism (Csv) is found to be most stable and 0.16 eV 
lower than the bicapped square antiprism (D4d). This ground state structure is consistent 
with previous LDA resultsil'ii. 

In Table IV, we compare the binding energy per atom Eh/n for Ge„ (n = 2 — 10) with the 
other theoretical and available experimental results. Due to the local density approximation, 
LDA calculation0ii has systematically overestimated the cluster binding energies. The more 
accurate binding energies for small germanium clusters up to five atoms has been provided 
by a more sophisticated G2(MP2) computationil. Although all the empirical parameters in 
our NTB model are fitted from dimer and bulk solid and there is no bond counting correction 
included, the experimental cohesive energies are fairly well reproduced by our calculation. 
Typical discrepancy between our calculation and experiment is less than 0.1 eV for those 
clusters. The successful description of binding energy within the present size range further 
demonstrates the transferability of the nonorthogonal tight-binding approach. In Table IV, 
we have also included the binding energies from Menon's NTB model for Ge„ clustersEl. 
Although the geometries of Ge„ found in their work is almost the same as our results, the 
binding energies of Ge„ starting from Ges iii Ref.l3 are about 0.10 ~ 0.2 eV higher than 
our results and experimental values. 

In summary, a nonorthogonal tight-binding model for germanium has been developed 
in this work. The transferability of the model is tested by various of bulk phases. The 
agreements between NTB model and ab initio results for bulk solids and small clusters are 
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satisfactory. For most Ge„ cluster with n = 3 to 10, the ground state geometries from tight- 
binding model coincide with those from ab initio calculation. The only exceptional cases are 
Ges and Geg, in which the ab initio metastable isomers are predicted as ground states by 
NTB scheme. However, the energy difference between the ground state configuration and 
the isomer is less than 0.01 eV/atom and within the accuracy of tight-binding approach. 
Therefore, the NTB model developed in this work can be applied to explore the configuration 
space of larger germanium clusters with n > 10, for which a global minimization at the 
ab initio level is significantly more expensive. Our further studies shall include a genetic 
algorithm for sampling the phase space and finding possible low energy structural isomers of 
germanium clusters. Thus, first principle structural optimization can be performed on these 
local minima structures. On the other hand, this model will be also employed to simulate 
the thermodynamic properties such as melting and growth process of germanium cluster, 
which require a long time scale in TBMD simulation. 

This work is partially supported by the U.S. Army Research Office (Grant DAAG55- 
98-1-0298) and the National Natural Science Foundation of China. The author (J.Zhao) 
are deeply grateful to Prof.E.Tosatti, Dr.J.Kohanoff, Dr.A.Buldum, and Prof. J. P. Lu for 
discussions. 
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